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We analyze a realistic microscopic model for electronic scattering based on the neutral differ- 
ential delay equations for point charges of the Wheeler-Feynman electrodynamics. We propose a 
microscopic model according to the electrodynamics of point charges, complex enough to describe 
the essential physics. Our microscopic model reaches a simple qualitative agreement with the ex- 
perimental results as regards interference in double-slit scattering and in electronic scattering by 
crystals. We discuss our model in the light of existing experimental results, including a qualitative 
disagreement found for the double-slit experiment. We discuss an approximate solution for the neu- 
tral differential delay equations of our model using piecewise-defined (discontinuous) velocities for 
all charges and piecewise-constant-velocities for the scattered charge. Our approximation predicts 
the De Broglie wavelength as an inverse function of the incoming velocity and in the correct order of 
magnitude. We explain the scattering by crystals in the light of the same simplified modeling with 
Einstein-local interactions. We include a discussion of the qualitative properties of the neutral dif- 
ferential delay equations and the boundary- value variational method of electrodynamics to stimulate 
future experimental tests on the possibility to complete quantum mechanics with electromagnetic 
models. 
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I. INTRODUCTION 



We discuss an attempt to complete quantum mechan- 
ics (QM) with an Einstein-local deterministic theory 
along with an analysis of realistic models for electronic 
diffraction. We introduce time-dependent interactions 
with the microscopic electronic trajectories inside the 
measuring apparatuses as a five-body-problem. We find 
that a microscopic model with delayed long-range inter- 
actions is enough to explain qualitatively the electronic 
scattering experiments. We analyze the microscopic 
models with the equations of motion for point charges 
of the Wheeler-Feynman electrodynamics^, whose com- 
plex qualitative behavior is worth stating clearly. The 
electromagnetic equations of motion are neutral differ- 
ential delay equations and have a qualitative behavior 
quite far from ordinary differential equations (ODE)^"— . 
We include a brief discussion of some generic qualitative 
properties of neutral differential delay equations relevant 
for the microscopic modeling. As with any modeling, 
if one starts with a narrow enough setup it is easy to 
obtain impossibility proofs, and the next step is usually 
harder, i.e., to define a wider minimum setup. We discuss 
our models in the light of available experimental work^ 
and a finite variational method for the Wheeler-Feynman 
electrodynamics^ . 

No single realistic dynamical model based on electro- 
dynamics has been analyzed down to microscopic detail 
up to now. The existing gap has two causes (i) electrody- 
namics itself was not a finished theory in the early 20th 
century; the first equation of motion for a point-charge^ 
appeared in 1938 and a form free of self-interactiortl only 



in 1945 and (ii) the theory of dynamical systems was not 
out in the early days of point-charge-electrodynamics, 
neither tools existed to understand the complex dynamics 
described by differential delay equations. On the other 
hand there is the popular quantum schism created by 
failed naive attempts to complete quantum mechanics 
with a Galilei-invariant mechanics based on ordinary dif- 
ferential equations with globally C°° trajectories and in- 
stantaneous interactions. Quantum mechanics (QM) was 
successfully developed as a Machian probabilistic theory, 
and there are several hints that a dynamics with New- 
tonian ordinary differential equations is not sufficient to 
complete QM. Even though the simplest quantum op- 
erators are pedagogically constructed by analogy with 
Galilei-invariant mechanics, a classical dynamics with 
ordinary differential equations is useless either to com- 
plete QM or to approximate the neutral differential delay 
equations of motion of point-charge-electrodynamics, as 
we show here. It is known that the no-go-theorems of Bell 
type fail for non-instantaneous interactions^ even with- 
out invoking the generic properties of neutral differential 
delay equations and variational methods with boundary 
conditions in the future, as discussed here. Nevertheless, 
there has been no attempt to complete QM with detailed 
modeling as regards the double-slit experiment. 

This paper is divided as follows: In Section |lT] we dis- 
cuss existing experimental results on double-slit scatter- 
ing, including a controlled ballistic experiment that found 
an unexpected qualitative disagreement^ and the later 
claims^ that the double-slit experiment should not dis- 
play the simple interference based on a single De Broglie 
wave. We give a simplified model based on electromag- 
netic trajectories with vanishing far-fields, a model that 
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explains the qualitative results of interference. In Section 
mil we discuss the equations of motion of point charges 
of the Wheeler-Feynman electrodynamics and the elec- 
tromagnetic variational method^. In Section Hill we also 
work out the details of a model for scattering with each 
slit represented by a heavy neutral atom, a system with 
little recoil that is not destroyed upon scattering, unlike a 
model for bowling. We introduce a simple approximation 
for the dynamics with (discontinuous) piecewise-defined- 
velocities and piecewise-constant-velocities for the scat- 
tered charge. Our model predicts the qualitative depen- 
dence of the De Broglie wavelength on the incoming elec- 
tronic velocity and calculates the De Broglie wavelength 
in the correct order of magnitude. In Section IIVI we put 
the discussions and conclusion. In appendix A we dis- 
cuss the qualitative properties of differential delay- and 
neutral differential delay equations with generic initial 
conditions. These qualitative properties when included 
in detailed modeling might offer further possible experi- 
mental tests for electromagnetic modeling and the com- 
pletion of QM. In appendix B we discuss some models 
for scattering by a periodic crystal. The modeling again 
reaches a qualitative agreement with the experimental 
results. 



II. ELECTRON SCATTERED BY A 
DOUBLE-SLIT 



Usually, the quantum explanation of double-slit 
diffraction is contrasted to simple billiard-ball-like- 
one-body-dynamics, and no microscopic model or 
time-dependent variables of the apparatus are ever 
introducediSiii. Let us start by defining this popular no- 
go modeling of scattering, henceforth called the billiard- 
ball-model; This model assumes the electron interacts 
with the slit walls by contact forces only. According to 
this model, if one closes one slit at a time, an outwardly 
scattered electron must have gone through the open slit 
without having interacted with the closed slit at all (be- 
cause interaction is by contact only)i^. Therefore, the 
statistical pattern at a far screen with both slits open 
should be the sum of two one-slit patterns, in bold con- 
tradiction with the predictions of Schroedinger's equation 
and the experimental results. On the other hand, the 
popular quantum model of double-slit diffraction is the 
scalar diffraction for the wave equation, a Schroedinger 
equation often pedagogically introduced by analogy with 
billiard-ball dynamics, with a potential vanishing every- 
where but at slit boundaries. We stress that even though 
it is time-proven that the analogy works well from the 
Galilei- invariant dynamics to abstract quantum opera- 
tors, using the analogy backwards to construct a dynam- 
ical system to complete quantum mechanics is far an- 
other problem. The reverse analogy from the scalar wave 
diffraction to the billiard-ball dynamics is not granted; 



The difficulties with this reverse analogy become appar- 
ent when we consider the realistic quantum mechanical 
model for a real electron with spin, approaching a realistic 
material slit composed of electrons and protons. 

The realistic quantum description is as follows; The 
incoming electron has a De Broglie distance of influence 
Xdb, so that its wavefunction overlaps the wavefunctions 
of the electrons at the material before it hits the slit. Be- 
cause of the exclusion principle, already at a distance 
XdBi the electronic spin could flip due to this overlap. 
Flipping the spin gives the incoming electron an orbital- 
angular-momentum-kick of the order of H (electronic spin 
carries an angular momentum of This change- 

by-spin-flip of the orbital-angular-momentum produces 
diffraction, even before the electrons hits any slit. Notice 
that quantum mechanics with spin suggests that the elec- 
tron suffers an angular-momentum kick without having 
to pass by any slit, actually even before going through one 
of the slits. Therefore, the careful quantum mechanical 
analysis of diffraction with spin suggests a microscopic 
model with an interaction that can change the orbital 
angular momentum at-a-distance (rather than a billiard- 
ball dynamics) . If the distance between slits is a and the 
incoming electron has mass m and velocity w, the largest 
orbital-angular-momentum about an axis passing by one 
slit belongs to a trajectory passing by the other slit, i.e., 
L = amv (higher trajectories simply hit the wall and 
do not pass to the other side). The largest scattering 
due to the h angular-momentum-kick of the flipping spin 
happens when L = amv ~ i.e., 

«=— , (1) 
mv 

which is the familiar slit-separation for significant scat- 
tering. In the electromagnetic microscopic modeling the 
electron interacts with both slits at a distance Xdb be- 
fore passing through one slit, unlike the interaction by 
contact of the billiard-ball model. So much to say the 
confusion generated by the billiard-ball straw-man-model 
is unacceptable moot. 

The first double-slit experiment^- was performed in 
1961, obtaining a scattering pattern in qualitative agree- 
ment with the quantum predictions. The first ballistic 
experiment with a single electron in the apparatus at 
a timei^ was performed in 1989 by Tonomura et ali^, 
obtaining a build-up of fringe patterns in qualitative 
agreement with quantum mechanics and the Huygens 
interference. It seems that so far no one has claimed 
quantitative precision beyond order of magnitude for the 
double-slit experiment. In 1994 an unexpected period- 
icity was found in a controlled ballistic experiment of 
double-slit scattering^. The anomalous periodicity was 
later explained^ as due to multiple scattering processes, 
but the quantum mechanical explanation made explicit 
use of a realistic model — . The explanation^ also made 
evident the need for a realistic quantum modeling and as 
a bonus seemed to rule out a precise experimental realiza- 
tion of the popular idealized experiment. This situation 
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should be contrasted to electronic diffraction in crystals 
and periodic structures, where Bragg directions can be 
determined very precisely. For the case of a crystal, the 
microscopic modeling captures the essential experimen- 
tal results even with ordinary differential equations, as 
discussed in appendix B. 

In the following we introduce a microscopic model for 
double-slit-diffraction consisting of the scattered electron 
in addition to an electron bound to each heavy positive 
charge at each slit site, which are separated by the slit dis- 
tance a. Our model uses the Wheeler- Fey nman equations 
of motion of point charges keeping only the far-field inter- 
actions. Our simplification stems from four reasons: (i) 
we disregard the short-range Coulomb interactions that 
are approximately cancelled by charge neutrality, (ii) we 
ignore the short-ranged Biot-Savart interactions that are 
small at low velocities, (iii) we keep the terms that in- 
troduce neutral delay in the differential equations, i.e., 
the far-field-couplings between point charges, which de- 
pend on the retarded/advanced accelerations, and (iv) 
the far-field interactions have the longest range, repre- 
senting naturally the first physical interaction suffered by 
the electron when approaching the double-slit and scat- 
tering away. We henceforth adopt a unit system where 
the speed of light is c = 1, the electronic charges are 
e = — 1 and we label quantities related to each elec- 
tron with indices k = 1,2,3. We define trajectories on 
a Lorentz four-space attached to an inertial frame 
by Einstein synchronization of clocks. A point in 
is defined by a time t and a spatial position x and we 
parametrize trajectories by time so that each trajectory 
is a function :x.k{tk) : K — > M^. As derived in Section Hill 
the low-velocity equation of motion of charge 1 in the 
far-field-only approximation is 

+^ X K3 X ^^-M)] 
X [nr3 X ^Mt^ )]- 

'13 

The equations of motion of the other charges are obtained 
analogously, by exchanging indices in Eq. ([2]). In Eq. 
@ unit vectors nf^^ = (xi(i) - Xalt^ )/|xi(t) - Xfc(tJ)| 
point respectively from the charge's advanced/retarded 
position Xfc(tfe) to the position of charge 1 at time t. In 
Eq. ^ the advanced/retarded times are defined by the 
implicit formulas 

tt=t±\Kk{t^)-^{t)l (3) 

and the respective distances in lightcone are defined by 

r±, = |xfc(t±)-x(t)|. (4) 



Notice that Eq. ([2]) is a neutral differential delay equation 
with advanced and retarded deviating arguments. In ap- 
pendix A we review some qualitative properties of delay- 
and neutral differential delay equations and discuss how, 
for generic initial data, solutions of delay equations must 
be defined piecewise. Moreover, our approximate equa- 
tions of motion depend linearly on the accelerations (see 
Eq. ©), so that arbitrary piecewise-constant- velocity so- 
lutions can be added to any balanced solution of the iso- 
lated system of bound charges. In Ref.— it is proved that 
if the bounded motion of two charges has vanishing far- 
fields, the dynamics of the bounded charges must have 
a sewing chain of velocity discontinuities. This gives rise 
to a natural forcing period defined by the sewing chains 
of the bound trajectories, as explained in the following. 

As discussed in appendix A, in general a neutral differ- 
ential delay equation like Eq. ([2]) has piecewise-defined 
solutions with discontinuous velocities and accelerations. 
A discontinuity in the acceleration of particle 1 is prop- 
agated in time by Eq. ([2]), as follows; In Figurel we 
illustrate a point po where the acceleration of charge 3 
is discontinuous. The equation of motion for charge 1, 
Eq. ([2]), involves the discontinuous past acceleration of 
charge 3 on the right-hand-side, yielding a discontinuous 
acceleration for charge 1 at point /3i(l). On the equa- 
tion of motion for charge 3 the discontinuous accelera- 
tion of charge 1 at point /3i(l) is on the right-hand-side 
and defines another discontinuity for acceleration 3 at 
point /3i(2). Successively, this produces a sewing chain 
of acceleration discontinuities at times /3i(fc), as illus- 
trated in Fig. 1, which define a stroboscopic sampling of 
each trajectory. The sewing chain involving trajectories 
1 and 2 is not illustrated, but it resonates and is part of 
the dynamics as well. A controlled ballistic experiment 
on double-slit interference found two wavelengths^; one 
near the De Broglie wavelength and another near half 
the De Broglie wavelength. Our electromagnetic model 
with state-dependent-delay naturally defines two wave- 
lengths by the four delayed interactions with the bound 
electrons, as follows; 

The equation of motion ^ and the respective equa- 
tions for the other charges have oscillatory solutions 
with decreasing wavelength when the electron approaches 
the inner gate of scattering. Low-velocity initial con- 
ditions approaching the inner gate from the left-hand- 
side with an exactly horizontal velocity have equal ad- 
vanced/retarded distances to the bound charges, i.e., 
rj'2 — r^2 — ^1*3 — '''13 = ^- I^or example, along tra- 
jectory 1 the bouncing times of the advanced interaction 
are times /3i(l), /3i(3), /3i(5) of Fig. 1, which in the 
far region are approximately separated by 2r and have 
decreasing time separations until the limiting period a 
for a perfectly symmetric horizontal velocity hitting the 
center of the inner gate (which never passes to the other 
side). In the near region the sewing chain of Fig. 1 starts 
to resonate also with the sewing chain between trajecto- 
ries 1 and 2 (not illustrated in Fig. 1), because they 
share the same frequency. The small vertical velocity 
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breaks the symmetry in the near region and the state- 
dependent distances to the bound charges start to differ, 
defining which sht the scattered charge passes through. 
The sewing chains for a shghtly non-horizontal trajectory 
3 excite an effective oscihation with period L superposed 
to the dynamics of the interacting bound charges. 

We henceforth assume that the scattered charge finds 
a resonant path along which the distances in lightcone 
to the bound charges are multiples of each other. The 
sewing chain of interactions with the bound charges ex- 
cite a period L smaller than the limiting value a of 
the symmetric centrally horizontal initial condition. In 
Section IIIII we estimate this resonant wavelength L by 
the distance of closest approximation between the scat- 
tered charge and the atom at the slit it passes through. 
The calculations assume the bound charges are bound to 
heavy positive charges sitting at each slit end, so that 
perturbations imposed by the incoming electron on the 
bound electrons are passed to the heavy charges at each 
site, producing smaller perturbations to disturb the scat- 
tered electron along its outer scattering path. 




FIG. 1: Illustrated is the trajectory of charge 3 (green) and 
trajectories of charges 1 (blue) and 2 (red), pictorially dis- 
played along the time direction, even though trajectories 1 
and 2 are bounded. The acceleration is discontinuous at point 
Po along trajectory 3, which is connected by the advanced 
interaction to point /3i(l) along trajectory 1, producing a 
forward sewing chain of breaking points, i.e., points /3i(2), 
/3i(3), /3i(4) and /3i(5). Illustrated is also a sewing chain of 
discontinuities starting from point po and proceeding in light- 
cone to trajectory 2, i.e., points /32(1), /32(2), /32(3), /32(4) 
and /32(5). 



The acceleration discontinuities along each bound tra- 
jectory oscillate periodically, so that an horizontally trav- 
eling third charge interacts with the discontinuity points 
along each bound orbit simultaneously (by symmetry), 
so that kicks received from each charge can cancel each 
other by synchronization. Therefore, the forward scat- 
tering direction is a solution by symmetry. Placing the 
asymptotically scattered constant- velocity in another di- 
rection on the right-hand-side of Fig. 2 is tricky because 
the scattered electron interacts with the discontinuities 



along each bound trajectory at different times. As dis- 
cussed in Section Hm the scattered velocity must jump to 
leave the Einstein-local momentum p7)) continuous. The 
options for scattering are the resonant Bragg directions 
illustrated in Fig. 2, where charge 3 interacts simultane- 
ously with the velocity discontinuities along both bound 
trajectories, so that condition (jI7p is automatically sat- 
isfied with a continuous constant velocity for charge 3. 
This simple microscopic model is in qualitative agree- 
ment with the experimental results on quantum interfer- 
ence. We stress that the sewing chain oscillations illus- 
trated in Fig. 1 are natural for our advanced/retarded 
equations, unlike the case of a model using an instanta- 
neous ODE. For example, the popular model using an 
ODE with instantaneous short-ranged Coulomb interac- 
tions is excluded, an in that respect even an ODE with 
long-range instantaneous interactions does not display 
oscillatory sewing chains. 




FIG. 2: Illustrated are the bound trajectories of electrons 
1 (solid blue) and 2 (solid red) at each material slit and the 
piecewise-straight-line trajectory of the scattered electron 3 
(solid green). Illustrated is the distance from particle 3 to 
particle 2 (green) and the distance from particle 3 to particle 
1 (blue). The different distances from electron 3 to the bound 
trajectories define a different delay for each interaction, which 
resonates when this difference is equal to the period L. 

Our simplified model explains the qualitative proper- 
ties of quantum double-slit scattering; namely (i) Ex- 
istence of favored lateral directions of scattering con- 
trolled by Bragg interference, which is due to the state- 
dependent delay, (ii) Existence of a forward direction of 
scattering by the double-slit system, and (iii) Scatter- 
ing with one slit closed at a time defines only a for- 
ward scattering peak while scattering with both slits 
open displays interference. The interference is due to 
the long-range nature of the forces and is present in any 
model with interaction-at-a-distance. The specific form 
of our interference turned out to be wavelike because we 
used the electrodynamics of point-charges, which inher- 
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ited state-dependent-delay from Maxwell's wave theory. 
Other models with long-range interactions should dis- 
play interference as well, but perhaps not precisely wave- 
like. In the next section we discuss a far-field model with 
piecewise-constant-velocities that predicts the inverse de- 
pendence of the De Broglie wavelength on the incoming 
particle velocity and yields a value in the correct order 
of magnitude. 



and 



B±(x,i) = T^x[- 



^(l±n±.v±)2^(l±n±.Vfe)3J' 

(9) 

where the advanced/retarded times are functions of t and 
X given by 



i± = i±|Xfe(t±)-x|, 



(10) 



III. PIECEWISE-DEFINED-SOLUTIONS AND 
FAR-FIELD-MODEL 

In the Wheeler-Feynman electrodynamics one does 
not solve differential equations for the electromagnetic 
fields, but rather the trajectories are the critical points 
(minimizers) of a variational method with mixed-type 
boundaries^. The electromagnetic fields are simply the 
coupling terms to the other particle's trajectories as de- 
fined by the Euler-Lagrange equations of the variational 
method. Rigorously speaking, these Euler-Lagrange- 
coupling-fields are defined only on the trajectories, even 
though given by the exact same usual formulas of 
Maxwell's electrodynamics^. Our model uses the equa- 
tions of motion of point charges keeping only the far- 
field interactions. We disregard the short-range Coulomb 
interactions that are approximately cancelled by charge 
neutrality and ignore the short-range Biot-Savart inter- 
actions that are small at low velocities. In a unit sys- 
tem where c = e = 1, the low- velocity Wheeler-Feynman 
equations of motion can be expressed in the familiar 



mfc^Xfe(t) 



-E(xfe,t) - Vfe X B(xfc,i), 



(5) 



where vj. = dy^\^{t)l dt. The electric far-field-coupling- 
term for each point-charge in the Wheeler-Feynman elec- 
trodynamics is given by a semi-sum of advanced and re- 
tarded fieldsi. 



£;(x,t) = ^E+(x,t) + iE-(x,i), (6) 

while the far-magnetic field is given by 

B(x,t) = in+ X E+(x,t) - x E-(x,t), (7) 

where unit vectors = (x — X2(i^)/|x — Xfc(t^)| point 
respectively from the charge's advanced/retarded posi- 
tion Xfc(ifc) to point X at time t ( see Chapter 14 of-**). 
The far-electric and far-magnetic field of each charge k 
are defined piecewise by the Lienard-Wiechert formulas^ 



E^(x,0 = 



(1 ± n 



(8) 



and the distances in lightcone are functions of t and x 
defined as 



Xfc(ife) -x| 



(11) 



where is defined by Eq. (ITU]) . In Eqs. ([5]) and 
([5]) = rfxfc/difeiji and = d^Xfe/dt^l^i are respec- 
tively the charge's velocity and acceleration at the ad- 
vanced/retarded time i^. It is possible to show that the 
derivative of the advanced/retarded time respect to t is 
given by 



dU_ 
~dt 



1 



(l±%±-v±)^ 



(12) 



using Eq. ([TU]) and the chain rule. Using Eqs. © and 
(jlOP with the chain rule we can re- write the far-magnetic 
field as 



B±(x,i) = T^x— x.(t±). 



while the far-electric field is given by 

E±(x,t)-±n±xB±(x,0. 



(13) 



(14) 



The far-field at each charge is the semi-sum of the ad- 
vanced and the retarded fields and of the other 
two charges. For example, the low-velocity equation of 
motion of charge 1, (Eq. ([S]) with /c = 1), is a simple- 
looking neutral differential delay equation 



d^ , , 
mi— xi(t) 



(n+ 



12 



X [n+ X ^X2(t+)] (15) 



("12 - ) 



(43 + Vp 



X ["12 X ^^^^2(^2 )] 
^2 



(n-|o — v-i ) . (P' , , , 

+ ^ " - X [7.13 X 

where the Xfe(t^) for fc = 2,3 are implicit functions of 
time t by Eqs. (fTO| with x = Xi(t). The equations of 
motion for the other (negative) charges are obtained by 
exchanging indices in Eq. p3|) . and Eq. © is the low- 
velocity-limit obtained by dropping the velocity terms of 
the right-hand-side of Eq. (fTS]). 

The fields of the heavy positive charges are essential 
to cancel the short-range Coulomb interaction and might 
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participate to create vanishing far fieldslS. In turn, the 
motion of the bound electrons, henceforth assumed pe- 
riodic, is essentially responsible to create a vanishing 
far-field that leaves the scattered electron undisturbed 
at large distances. For this, velocity discontinuities of 
the bound electrons are absolutely needed, as follows; 
As shown in^^, if a bounded motion of two charges pro- 
duces a net far-field that vanishes almost everywhere, 
then their trajectories must have discontinuous veloci- 
ties on sewing chainai^. This is the physical need for dis- 
continuous velocities. Vice-versa, the scattered electron 
can go out in a motion with piecewise-constant-velocity 
that produces no far-fields (because of Eqs. (|8]) and (O), 
leaving the bound charges undisturbed. An approxima- 
tion capturing the essentials of this five-body-dynamics 
is by orbits with discontinuous velocities'^ of the elec- 
tromagnetic variational method^. In the following we 
approximate the trajectory of scattered charge 3 by two 
segments with piecewise-constant-velocity, having a sin- 
gle point where velocity 3 is discontinuous. 

An orbit with a discontinuous velocity can be a mini- 
mizer of the variational method only if some other veloc- 
ity is discontinuous in the past/future lightcone, in a way 
that conserves an Einstein-local momentum^. For these 
minimizers, the Wheeler-Feynman equations of motion 
hold piecewise, and at velocity discontinuity points an ad- 
ditional condition for an extremum is the continuity of an 
Einstein- local-spin-like- four- momentunti^ (^i,Pi) with 



(^) 



-(c-n,. 



- (c + n. 



13+ • 



r 



(16) 



and vector component 



P, = 



(c -I- n,. 



(17) 



which must be continuous at each velocity discontinuity 
point of each continuous trajectory, despite the discon- 
tinuous velocities on both trajectories^^. Equations (fTT]) 
and (|16p are written in a unit system where the speed 
of light is c and the electronic charge is e, in order to 
make dimensional dependences explicit. In the following 
we use discontinuous trajectories as an intuitive approx- 
imation tool to estimate the dynamics without having 
to solve the detailed state-dependent neutral differential 
delay equations. 

In order to introduce a single discontinuity for the scat- 
tered velocity, our simplification further assumes that at 
its velocity-discontinuity-point charge 3 sees the single 



discontinuity of velocity 1 at a distance L and the discon- 
tinuity of velocity 2 also in lightcone at a larger distance 
nL. The velocity discontinuities along the bound orbits 
must cither see each other in lightcone (which requires 
a resonance), or we must introduce a larger number of 
velocity discontinuities. After the collision, the bound 
electron collides with the heavy positive charge, avoiding 
the destruction of the scattering system by moderating 
the size of velocity discontinuities. The remaining ex- 
cited oscillations have a subsequent influence on the out- 
scattered charge. In this situation we can produce the 
following estimate for the De Broglie wavelength; We use 
Eq. (|17|) to estimate the distance of closest approxima- 
tion L at the discontinuity point where charge 3 collides 
with the bound charge and its vertical velocity 3 jumps 
from zero to a value of the order of jvaj. The vertical 
component of Eq. (fTTl) yields the estimate 



'Valval 



L(c- n3i_ • vi_) 



(18) 



where we have assumed (c — n3i_-vi_) ~ (c-|-n3i+ • Vi+) 
and I Vi_ I ~ I Vi+ 1 . The collision is not necessarily a small 
perturbation on the bound electron's dynamics, and the 
denominator on the right- hand- side of Eq. ([TSl) can be 
large, as for example when vi_ performs fast solenoidal 
ping-pong motionsi^. Notice that Eq. ([T5| yields the 
experimentally tested De Broglie length inversely pro- 
portional to the incoming velocity. The bound electron 
subsequently collides with the heavy positive charge at 
the slit end (a proton of mass rup), sharing the momen- 
tum recoil and healing most of its velocity discontinuity. 
Otherwise the scattered electron destroys the scattering 
system and we are studying a bowling problem. We have 
at hand a five-body-problem with ten state-dependent de- 
lays of neutral type. Our approximation introduces a sin- 
gle velocity discontinuity along the scattered trajectory, 
which must be compensated by a velocity discontinuity 
in lightcone along the orbit of the bound electron. Veloc- 
ity discontinuities balance in lightcone by producing con- 
tinuous momenta (Eqs. (fT7| and (fT6|) ) at (the) breaking 
point(s) along the protonic trajectory and at (the) break- 
ing point (s) along the trajectory of the bound electron. 



I.e. 



eVc 



riC- Hpe- • Vi_ 



nieC 



r(c - He 



0, 
0, 



(19) 
(20) 



where r is the distance in lightcone at the subsequent 
collision and Vp is the (assumed small) protonic velocity, 
so that we ignore denominators depending on Vp in Eqs. 



in 

vi+) 



We further assume (c 
- (c-n3i_ -vi-) ~ c- 



■'■pe 



^ (c + n. 



■pe+ 



|vi_|, SO that eliminating 
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r from Eqs. p9|) yields 



_ .j^-^. 

(c- n3i_ • vi_) me 



(21) 



and the De Broglie length Xdb = L defined by Eq. p8|) 
is 



me c mslval 



(22) 



Formula ([22)) yields a length four times smaller than the 
popular double-slit formula, a rough estimate relying on 
an approximate solution with discontinuous velocities of 
variational electrodynamics. In our five-body model, the 
electron interacts with a disturbance it has created, much 
like the experiment of Ref.^^ 

Apart from a proportionality factor, the dependence 
h oc {mp/meY^^ singled out as the only power-law 
consistent with the hydrogenoid spectroscopic formula, 
as follows; The spectroscopic frequency of the hydrogen 
spectrum for the rii — )■ n2 transition is 



Wl2 



niee 



1 



1 



(23) 



If the protonic mass is multiplied by two, as for a deu- 
terium atom, our power 2/3 law predicts a four-times- 
larger H^, so that the spectroscopic frequency scales 
to another line of the hydrogen spectrum, i.e., the fre- 
quency of the 2ni — )■ 2n2 transition. Together with the 
estimates of the spectroscopic lines estimate (|22|) of 
the De Broglie length suggests that Wheeler- Fey nman 
electrodynamics generalized to discontinuous velocities 
can offer an alternative to complete quantum mechanics. 
The derivation of a fundamental formula like Eq. p2|) is 
foreign to the Machian mentality of quantum mechanics; 
Notice that even though Eq. was derived using a 

hydrogen atom at each slit end, another model using a 
heavier atom with n electrons and a proportionally larger 
nuclear mass nrup yields the same estimate because there 
are n electrons to share the recoil in Eqs. (|19p . A model 
using a large number of atoms along the inner gate of 
length a would also lead to the same estimate, so that 
formula dm is somewhat "universal". 



IV. DISCUSSIONS AND CONCLUSION 

The danger of replacing the Wheeler-Feynman equa- 
tions of motion^ by a Coulombian ODE, besides lack of 
Poincare invariance, are that these are only part of the 
conditions for a minimizer of the boundary-value varia- 
tional method with boundaries in past and future^^ . As 
for the fact that these two-body equations of motion seem 
to reduce "formally" to a Coulombian ODE at low ve- 
locities and large separations, the qualitative theory of 



appendix A teaches us that reduction to an ODE with 
C°° trajectories is not granted with a neutral differential 
delay equationi^. From the perspective of the variational 
method with boundaries in past and future-'^^, the set of 
boundary conditions leading to a many-body dynamics 
with C°° trajectories is probably a set of measure zero. 
As discussed in appendices A and B, this is certainly not 
the most general set of initial conditions and it is unlikely 
that initial conditions leading to C°° smooth motions 
with no velocity /acceleration discontinuities are ever go- 
ing to be met in experiments. More likely a larger class of 
initial conditions may lead to solutions where discontinu- 
ities do not explode asymptotically and stay bounded by 
laboratory walls ( i.e., respect the constraints of the mi- 
croscopic model). The fact that the electromagnetic vari- 
ational method has boundaries in past and future^'^^ is a 
warning for Bell-type proofs not to assume reduction to 
an ODE, at variance with the two-point-boundary-value 
problem of Hamilton's principle of instantaneous Newto- 
nian mechanics. Even reduction to piecewise orbits 
is an extra hypothesis that narrows down the qualitative 
possibilities of the boundary value variational method of 
electrodynamics^. We hope to enrich the attempts to 
complete QM and to prevent/enlighten attempts at im- 
possibility proofs based on "microscopic modeling with 
trajectories" . 

We have started on detailed modeling motivated by 
the subsequent works^ that tried to explain the unex- 
pected experimental result-. In the light of a detailed 
model with an interaction-at-a-distance rather than by 
contact, the contradictions of the billiard-ball-model dis- 
appear. Unfortunately, the subsequent explanations^ of 
the experimental disagreement^, seem to suggest that the 
double-slit experiment does not have a simple quanti- 
tative/qualitative precision, or it might have provided 
a testbed for microscopic modeling. Notice that the 
Wheeler-Feynman electrodynamics with neutral differen- 
tial delay equations is an Einstein-local theory of inter- 
actions; The velocity discontinuities along sewing chains 
involve only Einstein-local-two-body-events connected by 
a lightcone^, much like " Einstein- local-quantum-jumps" . 
The agreement of our far-field-and-piecewise-constant- 
velocities-model with the De Broglie wavelength is a 
test for the use of minimizers with discontinuous veloci- 
ties to approximate the solutions of the electromagnetic 
variational method*'. In fact, the success of our simple 
model of Section lllll to complete quantum mechanics sug- 
gests that the electromagnetic variational method^ is a 
physically sensible completion for the electrodynamics of 
point-charges. 

Several (if not all) of the popular failures to complete 
quantum mechanics with electrodynamics stem from the 
fact that electrodynamics of point-charges with C°° tra- 
jectories is an incomplete theoryi^. The popular reverse- 
analogy with quantum mechanics suggested in the 1930's 
an electrodynamics formally analogous to the instan- 
taneous Hamiltonian dynamics, a straw-man to gener- 
ate impossibility proofs and conundrums. The subse- 
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quent search for a quantization procedure starting from 
a Hamiltonian lead to the famous stall of the Wheeler- 
Feynman progran*2£ in the days before the no-interaction 
theorem^^. The no-interaction theorero^^i of 1963 was 
probably the first serious obstacle for the reverse-analogy 
approach in the days before the dynamical systems the- 
ory. 

The electromagnetic variational method^ with bound- 
ary data yielding globally orbits contains the reverse- 
analogy electrodynamics as the poor cousin. Ignoring 
the dynamical systems theory of appendix A, several 
formal embeddings of electrodynamics into a renormal- 
ized Hamiltonian dynamics used expansions of delayed 
arguments about non-generic C°° trajectories, expan- 
sions that are meaningless along orbits with discontin- 
uous velocities. Besides expansion of delayed arguments, 
Dirac's derivation of the self-force^ used Poynting's the- 
orem, which holds only for globally orbits^. The 
abundance of runaways of the Dirac electrodynamics of 
point charges again suggests that physics lies with the 
discontinuous-velocity-solutions of the boundary-value- 
variational-electrodynamics^. 

The leading motivation to construct the boundary- 
value-variational-method^ was to pave the way to study 
the electromagnetic-two-body-problem numerically. Ad- 
ditional tests for the variational electrodynamics with 
discontinuous velocities would be a microscopic model for 
Rutherford scattering. Also a detailed model for Comp- 
ton scattering on wax^* should include many electrons 
and a heavy nucleus, as discussed below Eq. (|23)) of Sec- 
tion lllll Further understanding on this neutrally-delayed 
dynamics should come after the development of a robust 
numerical integrator. A variational integrator shall be 
published elsewhere along with some numerical experi- 
ments. 



V. ACKNOWLEDGEMENTS 



Figure 1; Illustrated is the trajectory of charge 3 
(green) and trajectories of charges 1 (blue) and 2 
(red), pictorially displayed along the time direction, even 
though trajectories 1 and 2 are bounded. The accelera- 
tion is discontinuous at point po along trajectory 3, which 
is connected by the advanced interaction to point /3i(l) 
along trajectory 1, producing a forward sewing chain of 
breaking points, i.e., points /3i(2), /3i(3), /3i(4) and 
/3i(5). Illustrated is also a sewing chain of discontinu- 
ities starting from point po and proceeding in lightcone 
to trajectory 2, i.e., points /32(1), /32(2), /32(3), /32(4) 
and /32(5). 

Figure 2; Illustrated are the bound trajectories of elec- 
trons 1 (solid blue) and 2 (solid red) at each material slit 
and the piecewise-straight-line trajectory of the scattered 
electron 3 (solid green). Illustrated is the distance from 
particle 3 to particle 2 (green) and the distance from 
particle 3 to particle 1 (blue). The different distances 
from electron 3 to the bound trajectories define a differ- 
ent delay for each interaction, which resonates when this 
difference is equal to the period L. 

Figure 3; Illustrated in red is a initial condition for a 
scalar delay-differential-equation and the orbit (in black) 
after t = 0. The derivative from the left at the initial 
time i = is given by the history and in general differs 
from the derivative from the right along the continuous 
orbit. At t = 1 the solution is already . 

Figure 4; Illustrated in red is a initial condition for 
a scalar neutral-delay-differential-equation and the orbit 
(in black) after t — 0. The derivative from the left at 
t = is given by the history and in general differs from 
the derivative from the right along the continuous orbit. 
At t = 1 the solution still has a kink in the derivative, a 
discontinuity that generically propagates to the breaking 
point at t = 2 and to all successive breaking points. 
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VI. FIGURE CAPTIONS 



VII. APPENDIX A: DELAY- AND 
NEUTRAL-DIFFERENTIAL DELAY EQUATIONS 

Since the electromagnetic equations of motion are 
state-dependent neutral differential delay equations^, in 
the following we discuss briefly the qualitative properties 
of delay- and neutral differential delay equations^Ti^. We 
start with constant delay to simplify the exposition. A 
delay equation for a scalar function is defined as 

^=F{yit),y{t-l)), (24) 

where F{a, b) is usually a given C°° function. The right- 
hand-side of Eq. ([M)) depends of the present value y{t) 
and the past value y{t—l)^ so that to construct a unique 
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solution to Eq. (p4)) starting from t = one needs to 
provide y{t) on the entire segment t e (—1,0). In Fig- 
ure 3 a typical initial condition is illustrated in red 
while the ensuing continuous trajectory is illustrated in 
black. Notice in Fig. 3 that the derivative from the left 
at t = is given by the derivative of the arbitrary his- 
tory (red) while the derivative from the right at t = 
is F(y(0),y{—1)). The necessary initial history segment 
and a discontinuous derivative at the starting point are 
not seen in the simpler integration of ordinary differen- 
tial equations. For an ODE the initial value y{0) deter- 
mines a unique solution up to a maximum time tp where 
F{y) ceases to be Lipshitz-continuous. In solving Eq. 
([24l) with an arbitrary history there is a discontinu- 
ous derivative at the starting point, so that one can not 
impose Eq. (|24p at t — 0. Assuming the initial condi- 
tion is such that a solution exists, things are no so bad 
because Eq. (IM)) holds onwards, and moreover the kink 
at < = smooths out, as follows. At t = 1 when Eq. ([M)) 
first accesses < = as a past point, the left- and right- 
derivatives predicted by Eq. (|24p are equal precisely be- 
cause the trajectory is continuous, so that the solution 
y{t) is at t — 1. Successively, when the smoothed 
point t = 1 is first accessed by Eq. (l24t at t = 2 the 
derivative of Eq. ([M)) exists, so that the solution became 
C^, and so on, acquiring more derivatives at successive 
points where t is a multiple of the delay, henceforth called 
"breaking points"—""—. It is usual in solving delay equa- 
tions numerically- to restrict to derivatives from the left 
and derivatives from the right at breaking points. Last, 
for equations with state-dependent delay the breaking 
points depend on the solution, instead of being simple 
equally-spaced multiples of the constant delay. If the de- 
viating argument is a monotonically increasing function 
of time, the above smoothing is also the generic behavior 
expected for state-dependent-delay. 



y'(t)= F( y(t), y(t-1)) 




t=1 t=2 t 



FIG. 3: Illustrated in red is a initial condition for a scalar 
delay-difTerential-equation and the orbit (in black) after t — 0. 
The derivative from the left at the initial time f = is given 
by the history and in general differs from the derivative from 
the right along the continuous orbit. At t = 1 the solution is 
already . 



The integration of a neutral differential delay equation 
for a scalar function with generic initial data is different, 
as follows; We start again with constant delay to simplify 
the exposition^"^. A neutral differential delay equation 
for a scalar function is defined as 

^ = Giyit),yit~l),yit-l)), (25) 

where G{a, b, c) is a given function. The right-hand- 
side of Eq. (1^ depends on the present value y{t), on the 
past value y(t — 1) and on the past value of the deriva- 
tive y{t — I). A differential delay equation is said to be 
a neutral differential delay equation when the derivative 
depends on a past value of the derivative itself, as seen 
in Eq. (P5|). To construct a unique solution to Eq. ([25)1 
starting from t ~ again we need to provide y{t) for 
t e (—1,0) and moreover a initial history is needed 
because Eq. ([25)) uses the past derivative. As we already 
know, the constructed derivative is discontinuous, as il- 
lustrated in Fig. 4 for a typical initial history (red). 
In Fig. 4 the ensuing continuous trajectory is illustrated 
in black. Notice that when point t = is first accessed 
at t = 1 along the continuous solution, the derivatives 
predicted by Eq. ([25]) from the left and from right of 
t = 1 are different. We must therefore restrict to left and 
right derivatives at breaking points, so that solutions are 
defined only between breaking points. At t = 2 when 
point t = 1 is first accessed from the right along the 
continuous solution, the derivative discontinuity is still 
seen by G{y{t),y{t — l),y{t — 1)), unlike the case of Eq. 

In general for a neutral differential delay equation 
the derivative could stay discontinuous at all breaking 
pointa^Ti^. Surprisingly, having derivatives undefined on 
a set of measure zero is irrelevant for a variational method 
that involves an integral over the velocities^, but never- 
theless an important qualitative difference. So much for 
continuous solutions of one-component delay equations 
like Eqs. ([M]) and (^5]) . 

Again, the difference to neutral differential delay equa- 
tions with state-dependent delay is that breaking points 
depend on the solution, instead of being simple equally- 
spaced multiples of the constant delay. In the general 
case of neutral differential delay equation, there is the 
extra complexity of solution termination^ when the de- 
viating argument are not monotonically increasing func- 
tions of time. Otherwise continuation is possible with 
discontinuous derivatives. In the case of electrodynamics, 
deviating arguments are always increasingS, so termina- 
tion does not happen and continuation is always possible. 
In order to make the forward continuation well-posed, 
one needs a rule to calculate each velocity discontinuity, 
like for example our Einstein-local momentum continu- 
ity Eqs. (HZl) and (HH) of Section |llll The Wheeler- 
Feynman equations for the two-body problem are ex- 
pressible in the form (j25p using a 12-component vector 
y{t), particle positions taking up 6 components while par- 
ticle velocities take up the other 6 components. The 
smoothing properties of solutions of many-component 
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y'(t)= G( y(t), y(t-1), y(t-1)) 




FIG. 4; Illustrated in red is a initial condition for a scalar 
neutral differential delay equation and the orbit (in black) 
after t = 0. The derivative from the left at t = is given by 
the history and in general differs from the derivative from the 
right along the continuous orbit. At t = 1 the solution still 
has a kink in the derivative, which is general does not vanish 
and propagates to the next breaking point at t = 2 and to all 
successive breaking points. 



crystal through an instantaneous Hamiltonian interac- 
tion with a periodic potential. We consider a generic 
Hamiltonian ordinary differential equation (ODE) with 
a periodic potential, i.e., 



Vg exp(iG • x), 



(26) 



where summation is over reciprocal vectors of any generic 
periodic lattice. We use complex notation in Eq. ()26|) . 
with the condition that the potential is real, i.e., V-g = 
Vq, and e is included because we assume the potential 
is small. The dynamics of Hamiltonian (j26p is not in- 
tegrable in general but it can nevertheless be simplified 
with a canonical transformation designed to remove the 
potential, i.e., to transform the Hamiltonian into the free 
particle Hamiltonian to lowest order in e. The conditions 
for this simplification are as follows; We choose a sim- 
plifying quasi-identity canonical transformation^^, with 
generating function depending on the old coordinate x 
and new momentum P 

F{x,P)=x-P + e^FGeyip{iG-x), (27) 



differential delay equations are similar to the neutral dif- 
ferential delay case^S. For example, the 18-component 
retarded-only two-body equations of motion of Dirac's 
electrodynamics'' can have solutions along which some 
components are continuous (e.g. the positions) while the 
velocities and accelerations are discontinuous at all the 
subsequent breaking points^-^. 



VIII. APPENDIX B: ELECTRON (OR 
NEUTRON) SCATTERED BY A CRYSTAL 



We give an example that illustrates how a wave theory 
is not needed for the explanation of electronic diffrac- 
tion by crystals^, a simple but not so popular result. 
Unlike double-slit scattering, the interference found in 
crystalline scattering can be explained by a billiard-ball 
Hamiltonian model with instantaneous contact forces, 
even though certainly not a fundamental explanation. 
Next and in sake of an explanation compatible with the 
theory of relativity, we discuss modeling crystalline scat- 
tering with neutral differential delay equations, analo- 
gously to Section ini The simple model of Section HIl uses 
very little of electrodynamics and is based on a simple 
Einstein-local continuity of momentum. Any Einstein- 
local theory for nuclear interactions could be used with 
the same model to explain neutron diffractioniS. Several 
early works have explained diffraction by a crystal using 
the old quantum mechanics, which is based in Hamilto- 
nian mechanics^i2i. 

Our Hamiltonian model assumes that the incoming 
electron interacts with the bound electrons of a generic 



where the are so far arbitrary, respecting only the 
reality condition F-g = Fq- "^^^ old momentum p is 
related to the new P by 

_dF _ 
^ dx 



e'^^iGFc eyip{iG ■ x). 

G 



(28) 



Substitution of p given by Eq. (1^ into Eq. (^5)) yields, 
up to the first order in e, that the new Hamiltonian is 



H 



1, 



Vg + i{G ■ P)Fg\ exp(zG • x) + 0{(? 



(29) 

Hamiltonian II of Eq. ([^^ is still expressed in terms 
of the old coordinate x and the substitution x — X. 
into Eq. is only an O(e^) mistake. We want to fix 
the coefficients Fq to vanish the first-order terms of H 
as given by Eq. (j29p . so that the new momentum P is 
a constant of motion. First let us choose this constant 
vector P such that G • P 7^ for all vectors of the 
reciprocal lattice, a procedure yielding the Fq, by 



F, 



G 



-iVG 
G P' 



(30) 



and the new Hamiltonian, Eq. pQ]). becomes the free- 
particle Hamiltonian up to the first order in e, i.e.; 



H = \\P\^ + 0{e^), 



(31) 



having straight-line trajectories up to times of order 1/e^ 
and a constant P as already mentioned. The conclusion 
is that the trajectory is an approximate straight-line for 
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most initial conditions, i.e., if G • P ^ 0. These tra- 
jectories account for the forward beam, and constitute 
the most intense peak of scattering, in agreement with 
experiments. On the other hand, if the initial P is per- 
pendicular to a set of G's of the reciprocal lattice, it is 
not possible to remove the corresponding first-order term 
ofEq. (pS)) . In this case the simplified Hamiltonian has a 
first-order potential to scatter the particle. In two spatial 
dimensions, the set of directions satisfying G • P = for 
a given P are coUinear and the resonant G's are parallel, 
otherwise the following simplification is only an approx- 
imation; We further assume the Vg decrease sufficiently 
fast with \G\, so that we can take the Go with the least 
modulus as the essential unremovable term of the poten- 
tial, which yields 



H=-\Pf 



2e\VG^cos{Go■x), 



(32) 



a simple pendulum Hamiltonian. The higher-order terms 
produce crossings of the hyperbolic manifolds of the pen- 
dulum Hamiltonian in the usual way^'^'^^. The equation 
of motion for the P of the normalized Hamiltonian (l32l) 



IS 



dP 

— = [2e|VGjsin(Go-a;)]Go, 
dt 

from which it follows that the momentum kick is 



AP = Go 



OO 

j 2e|VbJsin(Go -0;)^^, 



(33) 



(34) 



a vector parallel to the reciprocal- lattice vector Go- An 
upper bound for the size of the momentum-kick is the 
maximum momentum excursion along the separatris22i^ 



lAPI 



so that we can re-write Eq. (j34]) in the form 



(35) 



AP: 



(36) 



The interested reader should read the literatur o^^i^^ on 
chaotic scattering with a Hamiltonian ODE. The result- 
ing momentum kicks are sensitively dependent on initial 
conditions as regards magnitude but come in a prescribed 
direction, i.e., along a vector of the reciprocal lattice. In 
conclusion; along these resonant initial conditions, the 
particle suffers stochastic momentum-kicks in the direc- 
tion of a single reciprocal-lattice vector^^ as defined by 
Eq. (IM)) . This is precisely the von Laue condition^^ 
AP = hG, in full qualitative agreement with the Bragg- 
von Laue wavelike scattering of X-rays by crystals^. Our 
simple model lacks an explicit quantitative size of the 
angular-momentum kicks, and to calculate the integral 



(1341) we need the trajectories and the next-order terms 
that perturb the simple integrable Hamiltonian (15^ to 
a chaotic dynamics. We could adjust Vg„ to obtain an 
experimental angular-momentum-kick of the mechanism 
to agree with Eq. ([36]) . i.e., lo — ^4e|VG„/|Go|, an ad- 
justing beyond our simplified modeling. 

We discussed the above Hamiltonian model as a sur- 
prise not often stressed, but even though in qualita- 
tive agreement with experiments, modeling with an in- 
stantaneous Hamiltonian ODE is hindered in principle 
by the non-interaction theorem^^-; The non-interaction 
thcorcm^^ states that the only two-body motion describ- 
able by a Hamiltonian ODE with manifest covariance by 
the Poincare group is the free-particle motion^S,. Essen- 
tially, Hamiltonian dynamics is instantaneous, and rela- 
tivity introduces the finite speed of light, which demands 
delay-equations of motion with solutions defined piece- 
wise, as discussed in appendix A. The no-interaction 
theorem is not a surprise; delay equations need a his- 
tory, and only free-moving particles have a straight-line- 
past-history that can be reconstructed from point-like 
ODE data. Otherwise, for generic initial data (or bound- 
ary data) and nontrivial interactions, we can not expect 
an ODE to guess an infinite-dimensional history. Since 
free particles move with constant velocity until a col- 
lision, the non-interaction theorem suggests that if a 
Hamiltonian description is of any use, the physical prob- 
lem at hand must have particles moving with piecewise- 
constant-velocities, like the scattered charge in our model 
of Section IIIII 

In sake of a model compatible with the theory of rela- 
tivity, in the following we discuss scattering by a crystal 
using neutral differential delay equations and the same 
model of Section |lT] (essentially a non-Hamiltonian busi- 
ness). Our model places one bound electron at each site 
of an infinite crystal. The explanation follows in the 
same way of the double-slit modeling of Section HT] Again 
we assume the incoming charge superimposes a resonant 
perturbation of period L on the bound orbits at each 
crystalline site, in the same way explained in Sections 
im and IIIII The perturbations are again transferred to 
the heavy on-site charges to avoid a bowling effect and 
to heal most of the velocity discontinuities. There re- 
mains a small oscillation, so that the out-scattered elec- 
tron suffers delayed periodic kicks from these synchro- 
nized oscillations excited at the crystalline sites. Let u 
be the unit vector along the scattered velocity and A-u 
the change of this unitary direction upon scattering. The 
advanced/retarded delay time ±Afj- of interaction with 
the different crystalline sites due to the site-dependent 
(Arj) distances of retardation is^^ 



At J = Ail ■ Arj , 



(37) 



in a unit system where the speed of light is c = 1. The fa- 
vored directions of Bragg scattering are those along which 
Atj = nL for an integer n. Using Eq. ([57)) and the usual 
properties of reciprocal lattice vectors^^ that 



12 



G ■ Arj = 27rn, 

we find that the favored (resonant) Au must be along a 
reciprocal lattice vector G, i.e., 

Au = ^G. (38) 



We see that again a simple Einstein-local microscopic 
model comes in full qualitative agreement with the 
Bragg-von Laue wavelike scattering by crystals^^. The 
estimate for the De Broglie length L is the same given in 
Section Uni 
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